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Abstract 

Nereid, a satellite of Neptune, has a highly eccentric prograde orbit with a semi- 
^ ■ major axis larger than 200 Neptune radius and is classified as an irregular satellite. 

Although the capture origin of irregular satellites has been widely accepted, sev- 
eral previous studies suggest that Nereid was formed in the circumplanetary disk of 
Neptune and was ejected outward to the present location by Triton. Our time-series 
photometric observations confirm that the spin is stable and non-chaotic with a pe- 
riod of 11.5 hr as indicated by Grav et al. (2003). The optical colors of Nereid 
are indistinguishable from those of trans-Neptunian objects and Centaurs, especially 
those with neutral colors. We also find the consistency of Nereid's rotation with the 
size-rotation distribution of small outer bodies. It is more likely that Nereid origi- 
nates in an immigrant body captured from a heliocentric orbit which was 4-5 AU 
away from Neptune's orbit. 

Key words: Satellite — Nereid — Solar system: general 



1. INTRODUCTION 

Nil Nereid is the second largest Neptunian satellite of 170 ± 25 km in radius (Thomas 
et al. 1991). It has a prograde orbit with a semi-major axis of 5.5 x 10^ km and an eccentricity 
of 0.75 (Jacobson 2009). Because of its large and highly eccentric orbit, Nereid is categorized 
as an irregular satellite. Irregular satellites are generally believed to have been captured into 
the Hill sphere of the host planets from heliocentric orbits (e.g. Bailey 1971). However, as an 



Based on data collected at Subaru Telescope, which is operated by the National Astronomical Observatory 
of Japan (NAOJ). 
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exception, the history of Nereid is still unclear. Goldreich et al. (1989) suggest that Nereid 
could originally have been a regular satellite. In this model, it formed close to Neptune in the 
circumplanetary disk and then was transported outward due to perturbations by NI Triton. 
The orbit scattering could also change Nereid's orbit from circular to eccentric. 

Nereid is comparatively bright {V ~ 19.2 mag at opposition; Schaefer et al. 2008) for an 
irregular satellite. Numerous ground-based observations have been carried out for measuring 
the rotational properties. Some of them show large variations of the magnitude exceeding 
~1.3 mag with a rotation period of 8-24 hr (Schaefer and Schaefer 1988; Williams et al. 1991). 
Schaefer and Schaefer (2000) argue that Nereid has a nonperiodic brightness variation with a 
total amplitude of 1.83 mag. We note that Nereid was located in vicinity of the galactic plane 
when these earlier data were collected and the apparent large variations could have been due to 
the crowded field stars. In contrast, Buratti et al. (1997) and Brown and Webster (1998) report 
a small amplitude of less than 0.1 mag. Also, measurements by Voyager 2 over a 12-day interval 
found no brightness variation greater than 0.15 mag (Thomas et al. 1991). Grav et al. (2003) 
performed relative photometry of Nereid with 0.003-0.006 mag accuracy using the CTIO 4-m 
telescope and obtained lightcurves with a peak-to-peak amplitude of 0.029 ± 0.003 mag and a 
rotation period of 11.52 ± 0.14 hr. 

Different measurements so far report different levels of the rotational brightness varia- 
tion. Schaefer and Schaefer (2000) suggest chaotic variations of the rotation from year to year, 
like SVII Hyperion (Klavetter 1989; Wisdom et al. 1984). Dobrovolskis (1995) propose that if 
Nereid was originally a regular satellite and has an elongated body with larger than ~1% of 
deviation from sphere, Nereid may be in chaotic rotation for any spin period of longer than 
about 2 weeks due to the tidal effect. On the other hand, the long-term variability of rotation 
is possibly produced by precession of the spin axis caused by the gravitational torques from 
Neptune (Schaefer et al. 2008; Alexander et al. 2011). While the possibility of chaotic rotation 
was excluded by Grav et al. (2003), the rotation state remains uncertain. 

This paper presents time-series photometric data of Nereid collected at the 8.2-m Subaru 
Telescope. We determine the amplitude and period of the rotational brightness variation with 
high accuracy. They give strong indications of the rotation state of Nereid; whether stable or 
variable. In addition, the rotation period is available for investigating the links to other small 
bodies. We use the results to discuss the origin and dynamical evolution of Nereid in the final 
stage of planet formation. 

2. OBSERVATIONS AND MEASUREMENTS 

Observations were conducted on 2008 September 1, 2, and 29 UT using the Subaru 
Prime Focus Camera, Suprime-Cam (Miyazaki et al. 2002), mounted on the Subaru Telescope. 
The Suprime-Cam consists of ten 2048x4096 CCDs, which are arranged in a 5x2 pattern with 
interchip gaps of ~15". It covers a 34' x 27' field of view (FOV) with a pixel scale of 0''20. 
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Each image was obtained with a 240-sec exposure using the VR-hand filter with the center 
wavelength of 0.6 fim and the band width of 0.2 fim. 

One FOV around Neptune was taken for 1-2 hours each night (Table 1). Neptune was 
placed at an interchip gap to reduce its bright light. Nereid was ~7' away from Neptune and 
its solar phase angle shifted from 0.55° to 1.34° during our observations. The sky motion, 
2.5-4.0 arcsec hr~^, was slow enough to be approximately viewed as a point source within the 
exposure. The typical seeing was 0."70 in the first night, 0."65 in the second night, and 0."60 in 
the last night. 

The data reduction was performed using IRAF^ and SDFRED2 (Ouchi et al. 2004) with 
the following processes: overscan subtraction, fiat fielding, distortion correction, sky background 
subtraction, and relative photometry. We estimate the flux via aperture photometry using the 
APPHOT task of IRAF. The aperture radius is 1.4 times that of a typical FWHM on a night 
that is optimum to obtain high S/N data for a point source (Howell 1992). We select 9-10 
background point sources listed in the USNO-Bl.O catalog (Monet et al. 2003) as reference stars. 
Those are required to be brighter than 19 mag in i?-band with B — R of 0.90-1.04 mag. The 
color range corresponds to ±0.03 mag from V — R = 0.44 mag for Nereid (Schaefer and Schaefer 
2000) using the color correlation of main sequence stars (Drilling and Landolt 2000). The light 
curves are given by the relative flux of Nereid to their total flux at each shot. The photometric 
accuracy reaches 0.001-0.002 mag. The relative sensitivity among CCDs is corrected using 
the Sloan Digital Sky Survey (SDSS) photometric database (Abazajian et al. 2009). We have 
confirmed a linear correlation of magnitude between the Suprime-Cam yi?-band and SDSS 
r-band. 

Three fields at different airmass including the 4-7 Landolt standard stars with different 
V — R colors (Landolt 2009) were taken each night. The coefficients of atmospheric extinction 
at Nereid's color, namely V — R = 0.44 mag, are measured. The photometric zero point is 
estimated by V magnitude of the standard stars with airmass and color corrections. It allows 
VR-hand flux to be converted into V^-band magnitude. 

3. RESULTS 

The calibrated photometric measurements are presented in Table 2. The apparent mag- 
nitude of Nereid is derived from the relative flux to the reference stars and their total V 
magnitude. The brightness variation of a rotating body with such a small amplitude could be 
attributed to the changes in the cross-section area and/or albedo inhomogeneity. If Nereid is an 
ellipsoidal body with a constant albedo over the surface, it produces a double-peaked sinusoidal 
lightcurve. Indeed, Grav et al. (2003) show that a simple sinusoid model gives a good fit to 
Nereid's lightcurve. We assume that Nereid is covered by surface with a uniform albedo and 

^ IRAF is distributed by the National Optical Astronomy Observatories (NOAO). 
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its lightcurve is doubly periodic. 

For determination of the lightcurve, a variation of the observed magnitude (mobs) due to 
the heliocentric distance (r in AU) and geocentric distance (A in AU) must be corrected. We 
use standardized magnitude (m) derived from mobs — Slog ]^Q[r A/900]. The rotational periodicity 
is analyzed by a Fourier-transform method called the Lomb-Scargle periodogram (Lomb 1976; 
Scargle 1982). The lightcurve model is represented as a first-order Fourier series formulation. 
In addition, we approximated the opposition effect as linear increasing with approaching the 
opposition. 

The synthetic lightcurve with a period P to be fitted is 

A A 
m(t) = cos {uj(t — to)) + sin {uj{t — to)) + ka + mg, (1) 

where t, u, a, k, mg are the observed time, frequency {u = 2'jt/P), solar phase angle, slope 
of the phase curve, and basement magnitude, respectively. The peak-to-peak amplitude A is 



given by the coefficients Ac and As as A = JA"^ + A"^. The time epoch to is evaluated by 



c 



tan(2a;to) = =^iv TTT' 

where is the number of data points. 

We obtained a best-fit model with P = 5.75 ± 0.05 hr, A = 0.031 ± 0.001 mag, and k = 
0.138 ± 0.002 mag deg~^. The mean residual of the fitting is 0.0026 mag, which is comparable 
to the photometric errors. Figure 1 shows the lightcurve combined with all the data points 
corrected for the opposition effect and the best-fit model folded with the oscillation period P. 
The rotation period (2P = 11.50 hr ± 0.10 hr) and amplitude (0.031 ± 0.001 mag) well agree 
with those presented by Grav et al. (2003), namely 11.52 ± 0.14 hr and 0.029 ± 0.003 mag. This 
represents the constancy of Nereid's rotation state among August 2001, August 2002 (Grav et 
al. 2003), and September 2008 (this study). We found no evidence of the rotation variability 
over 7 years. The rotation period is much shorter than the lower bound to induce the spin-orbit 
resonance or chaotic rotation, about 2 weeks (Dobrovolskis 1995). Nereid is observed to be and 
likely to remain in a constant rotation state. 

Under the assumption of albedo homogeneity, the brightness fluctuation with 0.03-mag 
amplitude constrains the body shape. A triaxial ellipsoid with semi-axes a > b > c rotating 
about the c axis draws a lightcurve with amplitude of 



1/2 



A = 2.5 log 2n^-2 ■ 2n > (3) 

y a-^ cos"^ a + sm d ) 

where a = a/h, c = c/h, and 9 is an aspect angle between the line of sight and the spin axis 
(Lacerda and Luu 2003). Given a simple prolate spheroid [i.e. c= 1), a < 2.0 in ^ > 15°. 

For a prolate body, the period of precession derived from tidal torque is represented by 

3(S-1) P„.' 



where Ppre, -Porb, and Prot are periods of the precession, orbit, and rotation, respectively 
(Dobrovolskis 1995; Schaefer et al. 2008). It is aheady known that Nereid has e = 0.75, Porb 
= 360.13 days, and Prot = 0.48 days (11.5 hr). Schaefer et al. (2008) suggest that Nereid's 
precessional period is ~8 or ~16 years. Assuming Pprc = 16 years, the body must have a 
> 4.3. This elongated shape is hardly acceptable to reproduce a small lightcurve amplitude, 
supporting the assertion that Nereid has a non-chaotic rotation with no precession. 

It is widely known that atmosphereless bodies exhibit an exponential increase in re- 
flectance at tiny solar phase angle, called an opposition surge. Schaefer et al. (2008) deter- 
mined the phase curves of Nereid using multi-year-combined photometric data on 246 nights 
from 1998 to 2006 (see also Schaefer and Schaefer 2000; Schaefer and Tourtellotte 2001). They 
have a sharp brightening exceeding 0.3 mag deg"^ down from a ~ 0.5°, represented by the 
best-fit model for surface scattering (Hapke 2002) as 



On the other hand, Grav et al. (2003) reported the surge slope of 0.14 ± 0.08 mag deg~^ at a 
~ 0.4°, which disagrees with the measurements by Schaefer et al. (2008). 

Figure 2 shows the combined phase curves and their slopes of Nereid presented by this 
work and those of previous studies. Both the magnitude and surge slope of our data are 
consistent with the model curve of Eq.(5). In contrast, the surge given by Grav et al. (2003) 
is too shallow. The discrepancy is most likely due to the small range of phase angles spanned 
by their observations (0.35°-0.46°). There is no doubt that Nereid has a steep opposition 
surge, indicating that the coherent backscattering mechanism is dominant on Nereid's surface 
(Schaefer et al. 2008). 

4. DISCUSSION 

From the large, eccentric, and inclined orbit of Nereid, it is natural to consider its origin 
as a captured body from outside the Neptune system. Nesvorny et al. (2007) investigated the 
capture of irregular satellites in the planetesimal disk by three-body gravitational reactions 
during encounters between the giant planets. This model well reproduces the orbits of Nereid 
and most Neptunian irregular satellites. Also, the long-term dynamical stability of Nereid's 
orbit has been confirmed (Holman et al. 2004). It is perfectly possible that Nereid has been an 
irregular satellite since it started orbiting around Neptune. 

Suetsugu et al. (2011) found that long-lived temporary capture of planetesimals in the 
prograde direction occurs at a limited range of initial orbital elements. The capture rate has a 
peak at e ~ 3 and E c^O. Here, e is the initial orbital eccentricity scaled hj h = Ph/ a, where Ph 
is the Hill radius and a is the semi-major axis of the planet; E is the energy integral described 
as E = {(? + i^)/2 — 3b^/8 + 9/2, where i is the initial orbital inclination scaled by h and b is 



m = 19.655 - 2.51ogio {l + 0.54[1 + (1 - e-^)/^]/[2(l + z)'^]} , 
^ = tan[a/2]/0.0134. 



(5) 
(6) 
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the difference in the initial semi-major axis between the planet and a planetesimal scaled by 
i?H (Nakazawa et al. 1989). For Neptune, these conditions correspond to e ^ 0.08 and b ~ 3.8- 
4.6 AU in i ranging from to e (e and b are non-scaled e and b, respectively). The source of 
Neptunian irregular satellites with prograde orbits is likely to be a planetesimal population in 
the region between 25 AU and 35 AU if they were captured at the current location of Neptune. 

Although most planetesimals were lost from such a region, a portion of them possi- 
bly survives as TNOs and Centaurs. After dissipation of the disk gas, the giant planets un- 
derwent orbital migration due to angular momentum exchange with residual planetesimals 
(Hahn and Malhotra 1999; Tsiganis et al. 2005). This led to extensive outward transport of 
planetesimals from 25-35 AU to the Kuiper belt (Gomes 2003; Levison and Morbidelli 2003). 
Levison and Morbidelli (2003) indicate that small bodies in the entire Kuiper belt were formed 
within ~35 AU and were transported outward as a result of Neptune's migration. Combining 
those studies, prograde irregular satellites of Neptune may have the same origin as TNOs and 
Centaurs. If Nereid derives from a captured planetesimal, this implication applies to it as well. 

It is worth comparing Nereid with TNOs and Centaurs in term of some physical proper- 
ties to examine this hypothesis. Nereid has optical colors of B — V = 0.71 ± 0.04 mag, V — R 
= 0.44 ± 0.03 mag, and V" - / = 0.72 ± 0.05 mag (Schaefer and Schaefer 2000). The colors of 
351 TNOs/Centaurs are provided by the MBOSS color database^ (Hainaut and Delsanti 2002). 
We use only those with uncertainties less than 0.05 mag for a comparison with Nereid. Figure 
3 and 4 show the color-color diagrams, i? — V" vs. V — R and B — V vs. V — I, respectively. 
One can see that Nereid is indistinguishable from TNOs/Centaurs in the color-color planes, 
supporting a common origin between them. It is notable that Nereid's colors are similar to 
those of neutrally colored bodies rather than red bodies. Brown et al. (2011) suggest that TNOs 
formed in the inner part of the primordial disk retain H2O, CO2 and poor hydrocarbon species 
on the surface which lead to a neutral color due to UV and high-energy particle irradiation. 
Nereid may have been originally formed at the region inside of ~20 AU and transported into 
near Neptune's orbit. 

We also focus on the distribution of rotation period as another indicator of similarity 
between Nereid and TNOs/Centaurs. The rotation properties of small bodies, excluding tidally 
locked satellites, are dominantly changed by mutual collisions (e.g. Harris 1979). The coUisional 
evolution of rotation period depends on the size distribution, impact velocity, and bulk density, 
which are unique among small-body populations. Therefore, small-body populations with the 
same source are expected to have a common rotation distribution. 

Interestingly, TNOs/Centaurs seem to have a characteristic pattern in the size-rotation 
distribution. Larger bodies rotate more slowly if they are greater than a few hundred kilometers 
in diameter (Sheppard et al. 2008). We evaluate the probability of this trend using a large 
lightcurve database of 102 TNOs/Centaurs compiled by Duffard et al. (2009). Figure 5 shows 

^ Minor Bodies in the Outer Solar System (http:/ /www.eso.org/"'ohainaut/MBOSS). 
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the distribution of their rotation periods against diameter derived by an assumed albedo of 
0.1. There appears to be a monotonia increasing trend beyond 200-300 km. We calculated 
the Spearman's rank correlation coefficient (Spearman 1904; Zar 1972) of the sample including 
33 objects with D > 250 km except for the tidally locked Pluto-Charon binary and Haumea, 
a rapidly rotating dwarf planet. The correlation coefficient is Vg = 0.53, meaning that large 
TNOs/Centaurs have a positive correlation between rotation and size at the significance level 
greater than 99.5 %. This trend is explained by the random accumulation effect of angular 
momentum through collisions (Harris 1979; Dobrovolskis and Burns 1984). 

The rotation period of Nereid is plotted on Figure 5. We found that Nereid perfectly 
agrees with the size-rotation correlation of TNOs/Centaurs. This agreement supports the 
idea that Nereid and TNO/Centaur populations share the same origin. The capture scenario 
needs no additional spin-change events to explain the present rapid rotation, unlike the other 
hypothesis based on the outward ejection from nearby Neptune. It is more likely that Nereid 
is derived from a captured body rather than a primordial satellite of Neptune. 

5. CONCLUSIONS 

We present highly accurate lightcurves of Nereid on two consecutive nights and one night 
after one month later in September 2008 at phase angles of 0.5°-1.3°. The periodic analysis 
shows that the rotation period of 11.5 ± 0.1 hr, peak-to-peak amplitude in brightness variation 
of 0.031 ± 0.001 mag. They completely agree with the measurements in Grav et al. (2003). 
On the other hand, the steepness of magnitude increases due to the opposition surge is 0.138 
± 0.002 mag deg~^ which is consistent with the phase curve model presented by Schaefer et al. 
(2008). The conclusions of this work can be summarized as follows: 

1. The period and amplitude of Nereid's rotation are constant in 2001-2002 and 2008. We 
found no evidence of variability of the rotation state. The rapid rotation with a period 
much shorter than 2 weeks also rejects the chaotic rotation state suggested by Dobrovolskis 
(1995). These facts indicate that Nereid has a stable rotation without a tidal despinning 
effect. 

2. Long-lived prograde captures of planetesimals by planets occur in a narrow range of the 
initial eccentricity and orbital energy (Suetsugu et al. 2011). If Nereid was a body captured 
from a heliocentric orbit, the source region is likely to be from ~25 AU to ~35 AU where 
it is suggested that TNOs have been formed (Gomes 2003; Levison and Morbidelli 2003). 

3. The optical colors of Nereid are consistent with those of TNOs/Centaurs. Also, Nereid 
accords closely with the size- rotation distribution of TNOs/Centaurs. The available ob- 
servations support the hypothesis that Nereid derived from a captured body and shares a 
common origin with the TNO/Centaur population. 

We thank David Jewitt for helpful comments. This study is based in part on data 
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Table 1. Time-series photometric observations for Nereid. 
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Fig. 1. Composite liglitcurve of Nereid with a period of 5.75 hr. The vertical axis shows normalized 
magnitude corrected for the opposition effect. The circles, triangles and squares represent the data obtained 
on 2008 September 1, 2, and 29, respectively. The dashed curve is the best-fit model with a peak-to-peak 
amplitude of 0.031 mag. 
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Table 2. Photometric measurements of Nereid. 
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Table 2. (Continued.) 



MJD* Airmass my^ (mag) 
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* Modified Julian date. 

t Apparent V-band magnitude. 
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Fig. 2. (a) y-band basement magnitude of Nereid's lightcurves standardized by heliocentrie/geocentric 
distances (mobs — 51og]^o[''^/900] ; see text) as a function of solar phase angle. The data are derived from 
this work (filled circles), Schaefer et al. (2008) (open triangles), and Grav ct al. (2003) (open squares). 
The dashed curve is the best-fit model of the opposition surge given in Schaefer et al. (2008). (b) Slopes 
of the measured phase curves of Nereid in mag deg^^. The symbols and dashed curve are the same as (a). 
The horizontal bar on each point shows the range of phase angle used in estimation of the slope. 
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Fig. 3. Color-color diagram of i? — F vs. V — R for Nereid (open diamond; Schaefer and Schaefer 2000) 
and TNOs/Ccntaurs (filled circles; Hainaut and Delsanti 2002). 
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Fig. 4. Color-color diagram oi B — V vs. V — I for Nereid (open diamond; Schaefer and Schaefer 2000) 
and TNOs/Centaurs (filled circles; Hainaut and Delsanti 2002). 
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Fig. 5. The size-rotation distribution of TNOs and Centaurs presented in Duffard et al. (2009). The 
dashed hue shows a hnear regression hne to the plot of logD {D is diameter) vs. rotation periods of those 
bodies with D > 250 km except for Pluto, Charon, and Haumea. An open diamond represents Nereid. 
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